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In the standard Big-Bang nucleosynthesis (BBN) model, the primordial 7 Li abundance is overes¬ 
timated by about a factor of 2-3 comparing to the astronomical observations, so called the pending 
cosmological lithium problem. The 7 Be(n,a) 4 He reaction, which may affect the 'Li abundance, was 
regarded as the secondary important reaction in destructing the 'Be nucleus in BBN. However, the 
thermonuclear rate of 7 Be(n,a) 4 He has not been well studied so far. This reaction rate was firstly 
estimated by Wagoner in 1969, which has been generally adopted in the current BBN simulations 
and the reaction rate library. This simple estimation involved only a direct-capture reaction mech¬ 
anism, but the resonant contribution should be also considered according to the later experimental 
results. In this work, we have revised this rate based on the indirect cross-section data available for 
the 4 He(a,n) 7 Be and 4 He(a,p)‘Li reactions, with the charge symmetry and detailed-balance prin¬ 
ciple. Our new result shows that the previous rate (acting as an upper limit) is overestimated by 
about a factor of ten. The BBN simulation shows that the present rate leads to a 1.2% increase in 
the final ' Li abundance compared to the result using the Wagoner rate, and hence the present rate 
even worsens the ‘ Li problem. By the present estimation, the role of ' Be(n,a) 4 He in destroying 7 Be 
is weakened from the secondary importance to the third, and the 'Be(d,p)2 4 He reaction becomes of 
secondary importance in destructing 7 Be. 

PACS numbers: 26.35.+C, 24.30.-v, 24.50. +g, 27.20.+n 


I. INTRODUCTION 

The discrepancy between the predicted primordial 'Li 
abundance and the astronomical observation persists as 
a fundamental pending problem in nuclear astrophysics 
studies jl]-|i]. As a powerful tool to study the early 
universe, the standard Big-Bang nucleosynthesis (BBN) 
model employs only one parameter r], the number-density 
ratio of baryons to photons. With the more accurate 77 
value determined from the astronomical observations 3 , 
the predicted primordial 7 Li abundance is still a factor 
of 2-3 higher than that observed in the galactic halo 
stars { 5 ]. It has been argued that such discrepancy may 
arise from the uncertainties in the thermonuclear rates 
for those reactions involved in BBN @,ll IQ. In the 
past two decades, great efforts have been devoted to re¬ 
duce these uncertainties. For example, Smith et al. [8} 
made a new evaluation of the reaction rates for the most 
important twelve reactions involved in BBN. Later on, 
Descouvemout et al. [9] re-analyzed ten key reactions by 
using the .ff-matrix theory. However, the 7 Li discrepancy 
still remains unsolved with these updated data together 
with the recent investigations |7l. I10l - (l2l | for those possible 
reactions affecting the 7 Li abundance. 

The final 7 Li abundance in BBN is contributed both 
from the directly synthesized 7 Li, as well as those from 
the 'Be EC decays and 'Be(n,p) 7 Li reaction, however, 


the relic 7 Li nuclei mainly come from the latter process 
because most of the directly synthesized 7 Li is destroyed 
through the 'Li(p,a) 4 He reaction immediately. The 7 Be 
production is determined by the balance between the re¬ 
actions which synthesize it and those destruct it. There¬ 
fore, it is essential to accurately determine thermonuclear 
rates for those reactions involving ' Be nuclide. The key 
synthesizing reaction for ' Be is the 3 He(a, 7 )' Be reaction, 
which has been studied very well by various experiments 
and theories [l3l4T6j . The reactions for destructing 7 Be 
are 7 Be(n,p) 7 Li and 'Be(n,a) 4 He, of primary and sec¬ 
ondary importance, respectively. For the 'Be(ro,p) 7 Li re¬ 
action, its cross section has been studied in detailed in 
a wide energy range from 0.025 eV up to 8 MeV @, U3], 
which covers entirely the BBN effective energy region. 
As for the 7 Be(n,a) 4 He reaction, it could play a non- 
negligible role in direct ' Be destruction [l8j]. However, up 
to now, the ' Be(n,a) 4 He reaction rate adopted in the cur¬ 
rent BBN simulations and the reaction rate library |l9j 
is still the very old Wagoner rate [20], which is a simple 
theoretical estimation involving only the direct-capture 
reaction mechanism, and also without information on the 
sources of data and error estimate. In this work, we have 
derived the thermonuclear rate of 7 Be(n,a) 4 He based on 
the available indirect experimental data of 4 He(a,n)'Be 
and 4 He(a,p) 7 Li reactions, with the well-know charge 
symmetry and detailed-balance principle [2l|. With this 
new rate, we have examined its impact on the primordial 
7 Li abundance with a BBN code. 
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II. DERIVATION OF 7 BE(n,«) 4 HE CROSS 
SECTION 


By so far, there is only one direct cross-section mea¬ 
surement for the 7 Be(n,a) 4 He reaction. The experiment 
was performed by Bassi et al. [22| in 1963 by using the 
reactor thermal neutrons. Based on this limited infor¬ 
mation and the theory of nonresonant reaction, Wagoner 
made the first estimation of this reaction rate, which will 
be discussed in detail in the next section. In this section, 
we will present the method to derive the cross section of 
' Be(n,a) 4 He by using the available indirect experimental 
data. 

About 30 years ago, King et al. [13 measured both 
cross sections of 4 He(a,?r) 7 Be and 4 He(a,p) 7 Li. Un¬ 
der the assumption of charge symmetry, i.e., the neu¬ 
tron and proton configurations in the compound 8 Be nu¬ 
clide to be identical, they calculated the total cross sec¬ 
tions of 4 He(a,?r) 7 Be based on their data measured for 
4 He(a,p) 7 Li via the following equation, 

jyno pn 1 

-L g -L g / \ 

°Vi = cr no + a ni = pp 0 a Po ± p Pl a Pl , (lj 

where a no and a Po are the cross sections leading to the 
ground states of 7 Be and 7 Li, respectively; and a ni and 
<j Pi are those leading to the corresponding first excited 
states. Here, Pt is the penetrability factor defined by [HI, 

M- 


Pt{E,R) 


kR 

F 2 {E,R) + G 2 (E : R) ’ 


( 2 ) 


where k is the wave number, R the channel radius, and F( 
and Gi the standard Coulomb functions. It was shown 
that the calculated a n agreed with their experimental 
data very well. 

Inspired by this idea, we have derived the cross 
section of 4 He(a,n) 7 Be based on those measured for 
4 He(a,p) 7 Li J23l . I25| . as well as those measured for 
'Li(p,a) 4 He [26| with the detailed-balance principle. The 
experimental cross-section data of 4 He(a,p) 7 Li and those 
derived from 7 Li(p,a) 4 He are listed in the first and sec¬ 
ond columns of Tabic HI 

In the astrophysical higlr-temperature environment the 
excited states of nuclei involved are thermally populated, 
which can also contribute to the total reaction rate [HI 
H9j. The first excited state in 7 Be is located at 429 keV, 
which is too high to make appreciable contribution in 
the total rate comparing to the ground state at BBN 
temperature region. Thus, it is appropriate to calculate 
the thermonuclear rate of 7 Be(n,ct) 4 He by taking only the 
ground state of 'Be nucleus into account. Relying on the 
first term of Eq. ©, we have directly derived the cross 
section of 4 He(a,n) 7 Be (g.s.) by utilizing those data of 
4 He(a,p) 7 Li (g.s.) listed in Table|l] It is worthy of noting 
that the ground-state contribution required is difficult to 
be extracted from the total cross section measured [23j] 


TABLE I: Experimental cross-section data collected for 
4 He(a,p) 7 Li (g.s.), and the corresponding ones derived for 
'Be(n,a) 4 He, in units of mb. The adopted uncertainties in 
energies (in units of MeV) of E a and E c . m . are ±100 keV and 
±50 keV, respectively l±I. f25i]. 


E a 

&(a,p) 

E 

&(n,a) 

Ref. 



0.0113 

8.4±8.5 

[27] 



0.0196 

10.7±10.7 

[27] 



0.0510 

13.2±13.2 

[27] 

38.23 

13.0±0.4 a 

0.124 

17.5±10.4 

[26] 

38.41 

14.9±1.4° 

0.214 

23.1±8.2 

[26] 

38.54 

24.2±2.0 

0.279 

39.0±11.1 

[25] 

38.96 

35.5±2.5 

0.489 

59.4±10.3 

[25] 

38.97 

29.9±1.5 

0.494 

50.0±8.2 

[23] 

39.44 

49.2±3.1 

0.729 

79.1±9.4 

[25] 

39.80 

59.9±3.0 

0.909 

91.6±8.3 

[23] 

39.94 

64.6±2.6 

0.979 

96.8±7.6 

[25] 

40.56 

30.5±2.5 

1.289 

41.5±3.9 

[2 5} 

40.99 

27.0±2.2 

1.504 

34.4±3.0 

[25] 

41.35 

23.4±1.2 

1.684 

28.2±1.6 

[23] 

41.61 

17.9±2.0 

1.814 

20.8±2.4 

125] 

41.95 

12.0±0.6 

1.984 

13.3±0.7 

T23] 

42.57 

6.5 ±0.3 

2.294 

6.6±0.3 

[23] 

43.04 

13.1±2.1 

2.529 

12.6±2.0 

[25] 

43.52 

12.0±0.6 

2.769 

11.0±0.6 

[23] 

44.32 

52.0±2.6 

3.169 

43.6±2.2 

[23] 

45.64 

36.5±1.8 

3.829 

27.1±1.4 

[23] 

46.67 

27.2±1.4 

4.344 

18.6±1.0 

[23] 

47.65 

22.7±1.1 

4.884 

14.3±0.7 

[23] 

49.49 

15.1±0.8 

5.754 

8.6±0.5 

[23] 


“The cross-section data are derived from those of 7 Li(p,a) 4 He 
with Eq. [Tj Since their incident energies have no errors quoted 
in the original paper [2(1 . and we assume the same errors as in 
Refs. [2g.l25|. 


for 4 He(a,n) 7 Be, and also these data were only measured 
down to £’ a =39.43 MeV which is much higher than those 
measured for 4 He(a,p) 7 Li. 

In Eq. ([2]) , the penetrability factor depends on the 
channel radius R , orbit angular momentum £ and in¬ 
cident center-of-mass energy E. In the treatment 
of King et al., a channel radius of R=A.l fm [i.e., 
-R=ro(A[/ 3 ±A.]/ 3 ) with ro=1.41 fm] was utilized in both 
n±‘Be and p± 7 Li systems. In the present calculation, 
the penetrability factor Pe(E,R) for p± 7 Li is calculated 
by an RCWFN code |30|, and that for neutron is calcu¬ 
lated using Eq. m with Fe and Gz values tubulated by 
Feshbach and Lax [3l|. Since the a particle is the spin¬ 
less boson, the wavefunction for two identical a parti¬ 
cles must be symmetric under interchange. However, the 
wavefunction for an £=odd state in 8 Be is antisymmetric 
by interchanging the two a particles. This implies that 
the compound state in 8 Be must have even parity for the 
incident a+a channel. We know that both ground states 
for 'Li and 7 Be have odd parity, and hence the relative 
orbit angular momentum t must be odd. Since the orbit 
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centrifugal barrier [ozl{i+\)\ increases with respect to l, 
the p- wave (A=l) capture will dominate both the n+'Be 
and p± 7 Li exit channels. For each energy points listed 
in TablelU we have calculated the 4 He(a,n) 7 Be (g.s.) re¬ 
action cross section, and the associated uncertainty is 
estimated by considering the uncertainties in both ro 
(in range of 1.1 1.5 fm |32[) and incident energy (±100 
keV 123 . l25j| ) . Finally, with the detailed-balance princi¬ 
ple [2l| the cross section of 'Be(n,a) 4 He has been de¬ 
rived as listed in the third and fourth columns of TablelU 
Here, three data points are listed for i£ c . m .<0.1 MeV, 
which are derived based on the experimental data |27] of 
7 Li(p,a) 4 He simply by the first term in Eq.|T| The associ¬ 
ated uncertainties are estimated by taking those of ro and 
incident energies into account. Alternatively, the low- 
energy data of 4 He(a,n) 7 Be have been simply estimated 
by linearly interpolating two data points of 4 He(a,p) 7 Li 
at 37.48 MeV [25| and 38.09 MeV [26}, and ultimately 
converted to those of ' Be(n,a) 4 He. And we find that 
the interpolated results agree well with these three data 
points at E c , m . <0.1 MeV within uncertainties. 

III. WAGONER’S ESTIMATION 

The thermonuclear rate of 'Be(n,a) 4 He has only been 
estimated by Wagoner j2Cj in 1969, which has been gen¬ 
erally adopted in the current BBN simulations and the 
reaction rate library [lj|. In Wagoner’s paper, the reac¬ 
tion rate was calculated by an analytical formula 

N a (av ) = 2.05 x 10 4 x (1 + 3760T 9 ). (3) 

For the non-resonant neutron-induced reaction, <av> 
can be expressed by [24j, 1H, [HI 

(av) = 5(0) ± 0.33125(0)T g 1/2 ± 0.064635(0)T 9 , (4) 

with 5 (0) the astrophysical S factor near zero energy. By 
comparing Eq.[3]and Eq.[U the values of 5(0) and and it’s 
second derivative 5(0) (wrt velocity v) can be obtained 
by assuming a negligible first derivative 5(0). Thus, the 
direct-capture cross section for neutron-induced reaction 
can be calculated by [24|, H}, 

_ 5(0) __ 5(0) ± \S(Q)E 


In this way, the cross section of 7 Be(n,a) 4 He is calculated 
as shown in Fig. [T| 

There is no any information on the sources of 
7 Be(n,a) 4 He cross-section data used in the Wagoner’s 
original paper [2C| . Before his estimation, there was only 
one measurement made by Bassi et al. [22] at thermal 
neutron energy. Here, we have figured out the Wagoner’s 
way of thinking. The total cross section of ' Be(n,a) 4 He 
was assumed to be composed of two contributions, i.e., 


ci+ 02 , with a i for the p -wave capture of ( n,a ) and 02 for 



FIG. 1: (Color online) Comparison of different cross sections 
of 7 Be(n,a) 4 He among present work and other two differ¬ 
ent origins (Wagoner’s work [2(J and TENDL-2014 evalua¬ 
tion [371). 


the s-wave capture of ( 71 , 7 a), respectively. The former 
obeys the law of (J 20 CU, and the latter obeys the law of 
uiocl/u [lH Hi} ■ According to this way, Wagoner made 
his estimation based on an upper limit of 07 < 0.1 mb and 
02=155 mb measured by Bassi et al. at thermal neutron 
energy. Actually, the p-wave dominates the cross sec¬ 
tion (or total rate) at the energy region of BBN relevant. 
Therefore, Wagoner actually gave us only an upper limit. 

The Wagoner cross sections and ours are compared 
in Fig. m It shows that the present results are over¬ 
all lower than those of Wagoner (upper limit). How¬ 
ever, our results are only derived based on the indi¬ 
rect experimental cross-section data (including both the 
direct- and resonant-capture contributions), further mea¬ 
surements [35|, H(| are proposed to acquire the direct ex¬ 
perimental data. In addition, the recent theoretical eval¬ 
uation of TENDL-2014 [37[ based on a TALYS calcula¬ 
tion is also compared, which is entirely different from the 
present results as shown in Fig. [T] 


IV. REVISED THERMONUCLEAR RATE 

The thermonuclear 'Be(n,a) 4 He rate as a function of 
temperature was calculated by numerical integration of 
our experimental cross sections using the EXP2RATE 
code by T. Rauscher [38]. Rate values, obtained as the 
arithmetic mean between the low and high limits asso¬ 
ciated with the uncertainties on both the cross-section 
data and incident energies, are given in Table [TTJ The 
present rate can be well parameterized (less than 0.4% 
error in 0.1-5 GK) by the following expression (e.g., in 
the standard format of Eq. (16) in Ref. H3): 
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N a (<jv) = exp(—17.8984 + 0.2711T- 1 - 23.8918T g ~ 1/3 + 62.2135T g 1/3 - 5.2888 T g + 0.3869T g /3 - 22.6197 lnT 9 ). (6) 


TABLE II: Thermonuclear reaction rates for 7 Be(n,a) 4 He in 
units of cm 3 s -1 mol -1 . The ratio between present rate and 
Wagoner rate is listed in the last column. 


T(GK) 

Present 

Wagoner 

Ratio 

0.1 

(9.6±8.3)xlO s 

7.7xlO t> 

0.13 

0.2 

(1.7±1.3)xl0 6 

1.5xl0 7 

0.11 

0.3 

(2.3±1.7)xl0 6 

2.3xl0 7 

0.10 

0.4 

(2.9±2.0)xl0 6 

3.1xl0 7 

0.09 

0.5 

(3.5±2.2)xl0 6 

3.9xl0 7 

0.09 

0.6 

(4.2±2.4)xl0 6 

4.6xl0 7 

0.09 

0.7 

(4.9±2.6)xl0 6 

5.4xl0 7 

0.09 

0.8 

(5.6±2.8)xl0 6 

6.2xl0 7 

0.09 

0.9 

(6.4±2.9)xl0 6 

6.9xl0 7 

0.09 

1.0 

(7.2±3.1)xl0 6 

7.7xl0 7 

0.09 

1.5 

(1.2±0.7)xl0 7 

1.2x10 s 

0.10 

2.0 

(1.7±0.4)xl0 7 

1.5x10 s 

0.11 

2.5 

(2.1±0.4)xl0 7 

1.9x10 s 

0.11 

3.0 

(2.5±0.5)xl0 7 

2.3x10 s 

0.11 

3.5 

(2.9±0.5)xl0 7 

2.7x10 s 

0.11 

4.0 

(3.2±0.5)xl0 7 

3.1x10 s 

0.10 

4.5 

(3.4±0.5)xl0 7 

3.5x10 s 

0.10 

5.0 

(3.5±0.5)xl0 7 

3.9x10 s 

0.09 


Our new rate is about a factor of ten smaller than 
the Wagoner rate in BBN temperature range. As dis¬ 
cussed above, Wagoner just presented an upper limit for 
this rate. Therefore, we propose here that the cross sec¬ 
tion rri of 7 Be(?r,a) 4 He at the thermal neutron energy is 
about 0.01 mb, one order of magnitude smaller than the 
previous upper limit [22| . 

V. BBN SIMULATION 

We have investigated the impact of our new rates 
on the BBN predicted abundances of D, 3 He, 4 He and 
7 Li by using a recently developed code [4(|. This 
code can calculate the reaction flux for every specific 
reaction at arbitrary time point. In this work, the 
recent values for cosmological parameters and nuclear 
physics quantities, such as the baryon-to-photon ratio 
T7=(6.203±0T37) x 10 _1 ° 41] and the neutron life-time 
t=880.3 s [13, have been utilized in our model. The 
number of light neutrino families N„=2.9840±0.0082 de¬ 
termined by CERN LEP experiment [43} supports the 
standard model prediction of N v = 3, which is adopted 
in the present calculation. The reaction network has the 
same size as in Ref. [40] with nuclei A <16 from n to 
16 O, and the relevant reaction rates are adopted from 
the literature [§, 9. [3 HU, HU HU . 


Here, two simulations have been performed with the 
Wagoner rate and our new rate of the ' Be(n,o:) 4 He re¬ 
action studied, respectively. It shows that the predicted 
abundances of D, 3 He and 4 He do not change apprecia¬ 
bly, while that for 'Li increases 1.2% by adopting the 
new rate. Therefore, the present rate even worsens the 
7 Li problem. In order to clarify the reason, we have per¬ 
formed the reaction flux HI] calculations. Figure [2] shows 
the flow passing though arbitrary reaction involved in 
the BBN network in a timescale of about 10000 s (long 
enough for BBN) with the present rate. The flux passing 
through the 7 Be(n,a) 4 He channel is about 10 -12 mol/g 
(marked by a dashed red arrow), which is about a fac¬ 
tor of ten weaker than the result using the Wagoner rate. 
Based on the present estimation, the role of 'Be(n,a) 4 He 
in destroying 7 Be is weakened from the secondary impor¬ 
tance to the third, and thus the 7 Be(d,p)2 4 He reaction 
becomes the secondary important reaction in destructing 
' Be. In addition, our calculation shows that the cosmo¬ 
logical 7 Li problem could be solved (within the observa¬ 
tion uncertainties) provided only if that the thermonu¬ 
clear rate of 7 Be(n,a) 4 He is about 180 times larger than 
the Wagoner rate under the condition of without chang¬ 
ing the rates for those rest reactions. Therefore, it is 
unlikely to solve the 7 Li problem by a further study of 
this reaction. 


VI. CONCLUSION 

We have revised the thermonuclear rate of 
'Be(n,a) 4 He, which was regarded as the secondary 
important reaction in destructing the 7 Be nucleus in 
BBN, and ultimately tried to understand the cosmo¬ 
logical 7 Li problem better. The present work shows 
that the previous Wagoner’s result is overestimated by 
about a factor of ten. The BBN simulation shows that 
the adoption of the new rate can not yield appreciable 
change to the final 7 Li abundance, only about 1.2% 
enhancement, which even worsens the 7 Li problem. The 
resolution for this mysterious problem might resort to 
other mechanisms or new physics beyond the standard 
models (e.g., see Refs. [4|| [46)] )■ Another possibility is 
that the current observational data might not exactly 
represent the primordial 7 Li abundance. The detailed 
discussion about how to solve this pending problem is 
beyond the scope of this paper. 
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